Metallic ferromagnetism without exchange splitting 
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In the band theory of ferromagnetism there is a relative shift in the position of majority and minority 
spin bands due to the self-consistent field due to opposite spin electrons. In the simplest realization, 
the Stoner model, the majority and minority spin bands are rigidly shifted with respect to each other. 
Here we consider models at the opposite extreme, where there is no overall shift of the energy bands. 
Instead, upon spin polarization one of the bands broadens relative to the other. Ferromagnetism 
is driven by the resulting gain in kinetic energy. A signature of this class of mechanisms is that 
a transfer of spectral weight in optical absorption from high to low frequencies occurs upon spin 
polarization. We show that such models arise from generalized tight binding models that include 
off-diagonal matrix elements of the Coulomb interaction. For certain parameter ranges it is also 
found that reentrant ferromagnetism occurs. We examine properties of these models at zero and 
finite temperatures, and discuss their possible relevance to real materials. 



I. INTRODUCTION 



In the Stoner model for metallic ferromagnetism jjj], a 
rigid shift in the positions of the majority and minority 
spin bands occurs. The difference in band energy between 
minority and majority spin electrons is: 



- (k] = A = Urn 



(1) 



Here A is the exchange splitting, m the magnetiza- 
tion and U the exchange interaction. In deriving the 
Stoner model from a tight binding model, U corresponds 
to an on-site interaction, because the shift in Eq.(l) is 
momentum- independent. The Stoner model has been 
widely used for the description of metallic ferromagnets 

More elaborate approaches using spin density func- 
tional theory derive the band structure of majority and 
minority spins taking into account electron-electron in- 
teractions within the local spin-density approximation 
The result is exchange-split bands where the ex- 
change splitting is in general not constant across the 
Brillouin Zone. Nevertheless, here again there is always 
an overall shift in energy of majority and minority spin 
bands with respect to each other. 

In this paper we consider models where there is no 
overall shift in the relative position of the energy bands. 
Can ferromagnetism still occur? The answer is yes, if 
there is instead a change in the relative bandwidth of 
majority and minority spin electrons upon spin polariza- 
tion. Figure 1 contrasts the situation in the Stoner model 
and in the models considered in this paper. The gain of 
energy in the ferromagnetic state arises from an effective 
mass reduction, or equivalently a bandwidth expansion, 
of the majority spin electrons, leading to a lowering of 
kinetic energy. 



The possible role of effective mass reduction in driv- 
ing the transition to ferromagnetism was first proposed 
in connection with the 'double exchange' mechanism to 
describe manganese oxides More generally, we have 
suggested in the context of describing metallic ferromag- 
netism with a single band generalized tight binding model 
HU , that effective mass reduction upon spin polarization 
may be an essential feature of all itinerant ferromagnets. 
It was pointed out @] that the anomalous drop in resis- 
tivity and the negative magnetoresistance that is usually 
found in ferromagnets could be explained as arising from 
a change in effective mass rather than a change in scat- 
tering time as usually assumed H. Clear experimental 
evidence for this phenomenon has been found in opti- 
cal properties of the colossal magnetoresistance mangan- 
ites H as well as of some rare-earth hexaborides Q and 
magnetic semiconductors ]lO| ]. Whether this is a univer- 
sal phenomenon that occurs in all metallic ferromagnets 
remains an open question. 

In the model that we considered originally ferro- 
magnetism arises from the combined effect of exchange 
energy and band broadening as the system orders. Such 
may well be the situation applicable to real materials. 
However, to better understand the various theoretical 
possibilities it is of interest to ask whether ferromag- 
netism could be driven just by band broadening in the 
absence of exchange energy. As we show in this paper, 
this is indeed possible, if the bandwidth for one spin ex- 
pands while the one for the other spin contracts. In con- 
trast, in the model considered originally where exchange 
energy was also present, both minority and majority spin 
bands could expand upon spin polarization. Still, both 
in the original model and in the ones considered here a 
signature of this physics is that an overall shift in optical 
spectral weight from high to low frequencies takes place 
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as spin polarization develops. Application of this model 
to describe the shift in optical spectral weight observed 
in EuBq is discussed elsewhere fyj . 



II. THE MODELS 

We consider a single band of width D, with density of 
states g(e) symmetric around the origin: 



,£>/2 

/ deg(e)e = 

J-D/2 



(2) 



The magnetization m and number of electrons n per site 
are given by 

n, I° /2 d eff (e)[/(e T (e))-/( Ci (e))] (3a) 

-D/2 



n = , D ' 2 de 5 (e)[/(e T (e))+/( ei (e))] (3b) 

-D/2 

with / the Fermi function. In the class of models under 
consideration here the quasiparticle energies are given by 

e CT (e) = (1 - 2j x J CT - 2j a J_ <r )e - <jD[^m + h]-fi (4) 

Here, h = H/D, with H an external magnetic field (in 
units of energy) , and the bond charge I a is given by 



I* = I '" deg(e)(-^-)f(e a (e)) 

J-D/2 U l l 



(5) 



The "exchange" interaction k gives an overall rigid shift 
in the band energies, as in the Stoner model. The new 
features here arise from the interactions j% and j% that 
give rise to band narrowing that depends on the bond- 
charge occupation for each spin I a , which in turn will be 
a function of temperature and magnetization. 



III. RELATION TO TIGHT BINDING MODELS 

Consider a tight binding model with on-site, nearest- 
neighbor, exchange and pair hopping interactions 

H = - t 2_j ( c h c j° + h -c.) + Uy^riji 

<Cij > , <r i 



<%3> 



<ij > ,<j,(t' 



(6) 



The interactions arise from the various matrix elements 
of the Coulomb interaction between orbitals at nearest 



neighbor sites. A mean field decoupling leads to the 
quasiparticle energies 

n , , J V ; r , J J ' m (V + zJ) 
e a (k) = [1 -/,(_-_) - I_ a {- + -)]e k - (j- — 



t t 



t t 



with 



(7) 

(8a) 
(8b) 

the bond charge and magnetization parameters. The 
bandwidth and band energies are 



m =< c- T Ci| > - < c\^Cii > 



Ck 



D = 2zt 

= -<5> 



(9a) 
(9b) 



with z the number of nearest neighbors to a site and S 
a vector connecting a site to its nearest neighbors. The 



./'= — -— (10a) 

(10b) 
(10c) 

In connection with itinerant ferromagnetism we have con- 
sidered in the past the model with U and J only 0] , and 
also the one with U, J and J' = J jTJ). In terms of the 
notation in Eq. (4), the first case corresponds to 
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3i = 32 = J 
k = u + j 



(11a) 
(lib) 



with j — zJ/D, u = U / D. In this particularly simple 
case the hopping amplitude gets renormalized equally by 
bond charge density of equal and opposite spin. When 
one includes the pair hopping term J 1 , and in particular 
for J' — J, we have 



Ji = J 
32 = 2j 
k = u + j 



(12a) 
(12b) 
(12c) 



so that in this case more of the hopping renormalization is 
due to electrons of opposite spin. As seen in Ref. 15, this 
case favors more the existence of ferromagnetism near the 
top and bottom of the band compared to the first case. 

Note that according to Eq. (10) the band renormal- 
ization parameter j\ from the same-spin electrons should 
be expected to be smaller than that of opposite spin elec- 
trons, j2) due to the effects of nearest neighbor repulsion 
and of pair hopping. In this paper we will focus on that 
situation. 
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IV. GROUND STATE PROPERTIES 

We parametrize the interactions in Eq. (4) as 
h = aj 

h = i 



(13a) 



(13b) 



with < a < 1. The case studied previously with only 
exchange interaction J corresponds to a = 1, and the 
case with pair hopping J 1 — J corresponds to a = 1/2. 
If the system is not completely polarized the chemical 
potential is determined by the condition 



with 



and 



e T (e F t) = £i(£Fi) 



deg(e) 

-D/2 



n = n-f + ni 



to — n-f — nj. 



(14) 



(15) 



(16a) 



(16b) 



the total occupation and magnetization respectively. We 
assume for simplicity a flat density of states 



1 D D 



(17) 



We have seen earlier that the properties of these type of 
models do not depend strongly on energy variation of the 
density of states, unlike in the Stoner model. 

The occupations of up and down spin bands at T = 
are given by 



so that 



£f<t 1 
D 2 



efcf n + am — 1 
IT = 2 



and the bond charge densities are given by 
1 — (1 — n — am) 2 



so that the condition Eq. (14) is 

/ ■ • \ / 1 , 2 , ■ . > l - to 2 - (l - n) 2 
Oi - n){i - nf - in + ]2 ) ^ >- -- 

and for the parametrization Eq. (13) we obtain 



(18) 
(19) 

(20) 

jfe-1 
(21) 



3 = 



2(1 - k) 



(1 + a)(l - to 2 ) - (3a - 1)(1 - nf 



(22) 



The conditions for onset of fcrromagnctism and for full 
spin polarization result from setting to — and m = n 
respectively in Eq. (22). 

The resulting phase diagrams for various values of a 
are shown in Fig. 2. For all cases, partially polarized fer- 
romagnetic regions occur predominant around the 1/2- 
filled band. Note that for exactly half-filling no full po- 
larization ocurs for k < 1. Also, note that for small 
values of a, spin polarization is easier near the edges of 
the band, and for large values of a near the half-filled 
band. For the value a = 1/3 (not shown), the bound- 
ary for onset of spin polarization is independent of n, at 
j/(l-fc) = 1.5. 

However, the phase diagrams in Fig. 2 need to be 
supplemented by the condition that the coefficients of e 
in Eq. (4) need to remain positive, i.e. 



1 - 2 3l l a - 2j 2 7_ CT > 



(23) 



which may not be satisfied for some values of k. While 
this condition is always satisfied at the points where onset 
of fcrromagnctism occurs, it is not so for partial or full 
spin polarization. For k — 0, condition (23) is equivalent 
to the constraint 



n < 1 — to 



(24) 



for n < 1, and the corresponding symmetric one for 
n > 1. In particular, condition (24) implies that full 
polarization in the absence of exchange, i.e. k = 0, can 
only be achieved for n < 1/2 or n > 3/2. Partial polar- 
ization of magnitude to = cn, < c < 1 , can be achieved 
only in the range n < 1/(1 + c) and the corresponding 
symmetric one for n > 1. 

For non-zero exchange k, full polarization can be 
achieved for any n if the condition 



- a ~ 2 



(25) 



holds, which in particular is true if k > 0.5 for any a. 
(By contrast, the Stoner model requires k > 1 to give 
rise to ferromagnetism.) Otherwise, full polarization is 
only achieved in the range 



n < 



1/2 



1 



(26) 



l-a 



and the corresponding symmetric one above half-filling. 
Note that as a approaches unity, condition (26) will be 
satisfied for any filling for arbitrarily small exchange k. 

Note also that the case a = 1 is a very singular point. 
For exchange k = 0, there is no ferromagnetism in the 
model as no change in relative occupation of up and down 
spin bands can occur. On the other hand, for any k ^ 
the criteria for onset and full polarization can be satisfied 
for suitable j for any value of occupation n. 
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V. FINITE TEMPERATURES 

Eq. (5) for the bond charge yields 

h-h= r 2 &5(e)(^)[/(e T (e))-/(ei(e))] (27) 

J-D/2 L '/ Z 

and using 

M + hh = J -^(h + h) + ^=^('t - h) ( 28 ) 

we obtain on expanding the Fermi functions on the right- 
hand side of Eq. (27): 



h-h = {h-h)[h-h)D I deg(e)(—f(- 

J-D/2 u l l Ce 



D/2 



df 



in the absence of exchange, k — 0. Thus the equation 
that determines the critical temperature is found by can- 
celling (if — ijj on both sides of this equation: 



D/2 



df 



1 = (h-h)D I deg(e)(—f(-^) (30) 

J-D/2 D / 2 

More generally, in the presence of exchange we consider 
also Eq. (3a) for the magnetization, and obtain the set 
of equations 

h -h = (h - 3i){h - h)G2 + kmG 1 (31a) 



with 



m = (j2 - ji)(/ T - /|)Gi + kmGa (31b) 



Gi = D f D/2 deg(e)(-L) l (-?L) (32) 

J-D/2 D/2 dt a 

and the condition for the critical temperature is obtained 
by setting to zero the determinant of the coupled Eqs. 
(31), yielding 

[1 - (32 - ji)G 2 ][l - fcGo] - (32 - h\kG\ = (33) 



We focus here on ferromagnetism without exchange, so 
that the T c equation is given by Eq. (30), or 



l = (.72-Ji)G 2 



For infinite temperature G2 is zero, and as T decreases 
it becomes positive. At a critical temperature, that in- 
creases when the difference (j 2 — ji) increases, Eq. (34) 
will be satisfied. The situation is analogous to the usual 
Stoner model where the critical temperature is deter- 
mined by the equation 



1 = fcGn 



(35) 



At low temperatures we can approximate the Fermi func- 
tion derivative in Eq. (34) as a ^-function, yielding 



1 



Mj2-.i1) 

l-2(ji+j 2 )/ CT v D 



(36) 



with (e Fa /D) given by Eq. (19). As T -> 0, I a is given 
by Eq. (20) and the T c equation (36) reduces to the 
condition Eq. (22) for onset of ferromagnetism (m = 0) 
for the case k = 0. One can also obtain an approximate 
analytic form for T c in weak coupling by expanding G2 
in Eq. (34) to one higher order than used to obtain Eq. 
(36). This is in contrast to the case j% = j'2, where the 
lowest order equation analogous to Eq. (36), 



1 



l-2j(7 T +/ i ) 



(37) 



(29) directly yields T c at low temperatures on expanding I a 



= ro _ 71 ( k B T 2 
a 3 D 1 



2(ji+j 2 )/° 



(38) 



with 7° =/ CT (T = 0). 

In fact, the integral G2 (Eq. (32)) will not be mono- 
tonically decreasing as the temperature increases, as Go 
is. For example, for a half-filled band the Fermi func- 
tion derivative as T — » approaches a 5-function at zero 
energy, and the integral G2 will vanish due to the extra 
factors of e. As the temperature is increased from zero, 
G2 will increase until a critical temperature is reached 
where Eq. (34) is satisfied. At a higher temperature 
G2 will start decreasing until a second temperature is 
reached where Eq. (34) is satisfied. Thus we can ex- 
pect to find re-entrant ferromagnetism in this model for 
intermediate values of the effective coupling (J2 — ji). 

The magnetic susceptibility above T c is obtained by 
taking the derivative of the magnetization Eq. (3a) with 
respect to the external magnetic field. For the case 
ji 7^ 32 the magnetic field dependence of the bond charge 
needs to be taken into account. We obtain for \ = dm/dh 



X = 2 



G - (32 - 3i)(G o9 2 - Gl) 



(1 - fcG )[l - (32 - ji)G 2 )] - k(j 2 - n )G\ 



(39) 



with Gi given by Eq. (32). In particular, Go is the mag- 
netic susceptibility per spin in the absence of interactions, 
and Eq. (39) reduces to the usual RPA form if y\ = j 2 - 
(34) H More generally, the susceptibility Eq. (39) diverges as 
1/(T — T c ) as T approaches T c given by the solution of 
Eq. (33). Parametrization of the susceptibility as 



X(T) = 



3(T - T c ) 



(40) 



defines the effective moment p e f f , which is temperature- 
independent within the Curie- Weiss law, which is often 
seen experimentally. 
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VI. NUMERICAL RESULTS 

We solve the mean field equations numerically for vari- 
ous cases. Figure 3 shows the magnetization versus tem- 
perature for n = 0.5, a = (so that j\ — 0) and various 
values of j 2 = j in the absence of exchange (k = 0). In 
this case the phase diagram Fig. 2a shows that onset of 
magnetization at zero temperature occurs for j = 1.6 and 
full polarization for j = 2. However note that at finite 
temperature ferromagnetism also occurs for j < 1.6 as 
discussed in the previous section. For j = 1.6 the mag- 
netization approaches zero as T — > 0, and for j < 1.6 
reentrant ferromagnetism occurs. For j = 2 the sys- 
tem becomes fully polarized as T — > 0, and for j > 2 
the model becomes unphysical because the condition Eq. 
(23) fails to be satisfied as the temperature is lowered. 
Note that for these cases the shape of the magnetization 
curve is very different than what is obtained from the 
Stoner model or from molecular field theory. 

The effective bandwidth for spin a, D*, or equivalently 
the effective mass m*, are given by the relation 



D 

D* 



1 - 2 3l l a - 2j 2 /_ 



(41) 



with D and m the (spin-independent) bare bandwidth 
and mass respectively. Figure 4 shows the temperature 
dependence. For the majority spins, the effective mass in- 
creases as the temperature decreases, and then decreases 
as T is lowered below T c . For the minority spins, m* 
increases more sharply as T is lowered below T c . 

For small values of j however the effective mass for 
the majority spins also increases as the temperature is 
lowered further, and becomes eventually larger than its 
value at T c . One may ask why it is still advantageous 
to spin-polarize, when the energy will no longer be lower 
than what it was in the unpolarized state at T c . The 
answer is of course that at this lower temperature if the 
system was unpolarized the bandwidth for both up and 
down spin electrons would be much smaller, giving rise 
to a larger energy than in the polarized state (see Fig. 1 
of Rcf. 5(b) for an illustration of this effect). 

The behavior of the effective moment, Eq. (40), versus 
temperature for this case is shown in Fig. 5. For small j 
the temperature dependence is strong and it increases as 
T approaches T c , for larger j it decreases somewhat as T 
approaches T c . 

Results for a case where the band renormalization for 
the same spin, ji, is not zero, are shown in Fig. 6. The 
case a = 0.5 shown corresponds to parameters J' = J 
and V = in the tight binding model. Here the effec- 
tive moment is much less temperature dependent and no 
reentrant behavior is found for this set of parameters. 

Fig. 7 shows the effect of exchange on the behavior 
of the magnetization and effective moment. With in- 
creasing k the magnetization curve becomes steeper and 
resembles more the conventional behavior. The effective 
moment is rather constant with temperature for these 



cases corresponding to a Curie- Weiss law for the suscep- 
tibility. As discussed in Rcf. 5, for the pure Stoner model 
(exchange only) a strong temperature dependence of the 
effective moment is found. 



VII. OPTICAL ABSORPTION, 
MAGNETORESISTANCE AND 
PHOTOEMISSION 

The Drude formula for the real part of the optical con- 
ductivity is 



CTl(w) 



ne 



to* 1 + 0J 2 T 2 



(42) 



which describes optical absorption through intra-band 
processes. The conductivity sum rule that results from 
this formula is 



A ( \ 

aui<ji(Lu) = — 

v ; 2 to* 



(43) 



The cutoff ut m is to exclude transitions to other bands, 
which are processes not described by the expression Eq. 

(42) . We assume that iv m T » 1 in order to derive Eq. 

(43) from Eq. (42). Ordinarily, if no change in effec- 
tive mass occurs upon spin polarization, the optical ab- 
sorption given by these formulas will not depend on spin 
polarization (assuming a constant n). While one may ex- 
pect changes in optical absorption as function of temper- 
ature due to changes in the relaxation time t, no changes 
would be expected as function of magnetic field at fixed 
temperature. 

In contrast, in the models considered here we have 



and 



l+0J 2 T 2 TO* TO* 



Jo 2 y m* my 



(44) 



(45) 



and changes in optical absorption will occur both as a 
function of temperature and of magnetic field due to 
changes in the degree of spin polarization and corre- 
sponding changes in the effective masses. We assume that 
the effective mass is proportional to the inverse band- 
width, so that 



m 
m% 



= 1 - 2jJ a - 2j 2 I_ a 



(46) 



Figure 8 shows the typical behavior expected in low fre- 
quency optical absorption, for parameters j\ = 0, j 2 = 2, 
n = 0.5, k = as an example. The optical absorption 
decreases for T above T c as the temperature is lowered, 
and increases again rapidly below T c . When a magnetic 
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field is applied the optical absorption increases at all tem- 
peratures, with the largest increase at T c . Figure 9(a) 
shows the optical weight versus temperature illustrating 
this behavior. In figure 9(b) optical weights for majority 
and minority spin electrons are shown separately; when 
a magnetic field is applied, the optical absorption of ma- 
jority spins increases and that of minority spins decreases 
but less strongly, giving rise to the change seen in Fig. 
9(a). 

The dc conductivity is given by 



ne 2 T 



(47) 



and the magnetoresistance in our model will be deter- 
mined by the change in effective mass, or bandwidth, 
with spin polarization. Assuming a constant relaxation 
time the magnetoresistance is given by 



Ap p{H) - p(Q) 
P P(0) 

= n a (H)(l - 2j 1 J g (g) - 2j 2 J- g (ff)) 
£ ff n ff (0)(l - 2ji/„(0) - 2j 2 /_ CT (0)) 



-1 (48) 



and is shown in Fig. 10 for one case. It is maximum at 
T c and remains large well above T c . Similar results are 
found for other parameters in the model. 

Note that under the assumption that there is no sig- 
nificant change in the relaxation time with spin polar- 
ization, there should be a definite relation between the 
optical weight, 



W(T,H) = —j / dujax{uj) 



(49) 



and the magnetoresistance. For example, the quantity 

Ap/p(0) 



f(T,H) 



AW/W(H) 



(50) 



should be a constant independent of temperature and 
magnetic field. Deviations from constant behavior would 
indicate a dependence of the relaxation time on spin po- 
larization. 

Another important experiment that can shed light 
on the nature of the ferromagnetic transition is angle- 
resolved photoemission. For example, some results on 
manganites have recently been reported ]T^ | . In the mod- 
els discussed here the quasiparticle dispersion changes 
with magnetization, and hence one would expect char- 
acteristic signatures in the photoemission spectrum as 
function of temperature or magnetic field. Of partic- 
ular interest would be to obtain high resolution angle- 
arid spin-resolved photoemission spectra of ferromagnetic 
metals [|l7). Figure 11 shows an example of predicted 
angle-resolved photoemission spectrum for our model for 
one case without exchange: for majority spins, the Fermi 
velocity increases as the system polarizes and the quasi- 
particle peak disperses faster, and the opposite occurs for 



minority spins. For the non-spin-resolved spectrum (Fig. 
11c), two peaks develop as the system polarizes, one of 
which disperses slower and the other one faster than in 
the unpolarized case. However, for other parameters in 
the model, in the presence of exchange, one finds that 
both majority and minority spins disperse faster in the 
ferromagnetic state, reflecting the gain in kinetic energy 
upon spin polarization. 



VIII. SUMMARY AND DISCUSSION 

We have considered a model for ferromagnetism where 
the magnetic order arises from a modification of the 
width of the bands upon spin polarization instead of the 
usual exchange splitting. Thus, ferromagnetism arises 
from a gain in kinetic rather than potential energy. This 
represents the opposite extreme of the conventional un- 
derstanding of ferromagnetism as arising from a gain in 
potential energy (exchange energy) despite an associated 
cost in kinetic energy [[L8|. More generally, the models 
considered here also include an exchange splitting param- 
eter, in the presence of which ferromagnetism arises from 
a combination of gains in potential and kinetic energy. 

As we have seen, these models can be derived from 
a single band tight binding model that includes off- 
diagonal matrix elements of the Coulomb interaction. 
The terms that give rise to kinetic energy gain upon 
spin polarization originate in nearest neighbor matrix el- 
ements conventionally termed "exchange" and "pair hop- 
ping" . Physically, they represent " bond-charge repul- 
sion" i.e.' the repulsive Coulomb energy of elec- 
trons at the interstitial region between neighboring atoms 
rather than at the atoms themselves. In the present con- 
text the use of the term "exchange" or "Heisenberg ex- 
change" for such matrix elements is somewhat mislead- 
ing, because the driving force for ferromagnetism in fact 
is direct Coulomb repulsion of electrons in the bonds. 
The model discussed in this paper with exchange param- 
eter k = describes ferromagnetism without exchange 
splitting, even though the main interaction is the param- 
eter J which is conventionally called an exchange inte- 
gral. 

It is generally believed that band degeneracy is essen- 
tial to metallic ferromagnetism jL8|. It is certainly true 
that ferromagnetism usually arises is systems that have 
atoms with valence electrons in degenerate orbitals. Nev- 
ertheless, it is difficult to see how atomic orbital degen- 
eracy could be essential for example in the case of Ni 
and Ni-Cu alloys, which have a small fraction of a d- 
hole per atom and negligible polar fluctuations. Further- 
more, it should be kept in mind that all energy bands 
in metals are in fact non-degenerate, except at sets of 
states of measure zero in the Brillouin Zone (points or 
lines). Hence we believe that an approach that focuses 
on a single non-degenerate band, as in this paper, is sen- 
sible, whether that band arises from degenerate or non- 
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degenerate atomic orbitals. On transforming from Bloch 
to Wannier states for that band and computing the ma- 
trix elements of the Coulomb interaction with the Wan- 
nier states, the matrix elements discussed in Sect. Ill 
(as well as others) result. The question of what is the 
magnitude of these matrix elements in particular mate- 
rials is a difficult one. Within thepoint of view of this 
paper (and our previous work |^,^9|), band degeneracy 
may be important in determining the magnitude of these 
matrix elements, but not in determining the structure of 
the theory. 

The model considered here naturally gives rise to par- 
tial spin polarization, even with a constant density of 
states, unlike the Stoner model. In the absence of ex- 
change splitting, we found that ferromagnetism does oc- 
cur, however full spin polarization can only be achieved 
in a limited range of parameter space (far from the half- 
filled band). As the exchange splitting increases from 
zero, full polarization can be achieved over an increasing 
range of band filling. The condition for onset of spin po- 
larization may become more or less stringent as the band 
filling increases depending on the ratio of the same-spin 
to opposite-spin band renormalization parameters, 31/32- 

The magnetization versus temperature showed uncon- 
ventional behavior in the absence of exchange, with the 
curves being substantially less steep. An interesting fea- 
ture is that re-entrant ferromagnetism can occur in cer- 
tain parameter ranges. The magnetic susceptibility ver- 
sus temperature showed Curie- Weiss behavior for a wide 
range of parameters, although deviations can also occur. 

Re-entrant ferromagnetism, and in particular a situa- 
tion where the magnetic order increases as the temper- 
ature increases, is somewhat counterintuitive but by no 
means unphysical. In fact, experimental observations of 
such behavior have been reported in Y2N17 |^(| and in 
ThFe3 plj . Theoretical models that have been found 
to exhibit such behavior are Ising models with random 
ferromagnetic and antiferromagnetic exchange p2| , and 
even the Stoner model for special forms of the density of 
states However in the Stoner model such behavior 
is always associated with first order transitions. In the 
models discussed here, the re-entrant behavior can be 
easily understood. Consider a case where the magnetism 
has disappeared at sufficiently low temperatures. As the 
temperature is raised, the contribution of entropy to the 
free energy of the system increases relative to that of en- 
ergy. Now in our model (in the absence of exchange), 
by spin polarizing one of the bands will narrow and as 
a consequence the entropy contribution from electrons in 
that band will increase, so that this effect added to the 
energy gain from the electrons in the broadened band can 
lead to an overall decrease of the free energy upon spin 
polarization at higher temperatures. As the temperature 
increases further there will be a point where entropy from 
both bands will favor a vanishing of the magnetization as 
in the usual cases. 

The most characteristic feature of the models discussed 
here, however, both in the absence and in the presence of 



exchange splitting, arises in optical properties. A trans- 
fer of optical spectral weight from high frequencies to 
low frequencies always occurs upon spin polarization in 
the presence of the interactions j± and 32 that modify 
the bandwidth. Furthermore, above T c these interactions 
cause spectral weight to be transfered from low to high 
frequencies as the temperature is lowered. These phe- 
nomena have been observed in manganites || and hex- 
aborides Q as a function of temperature. Within the 
models considered here, the same effects should be seen 
as a function of magnetic field. In fact, recent optical ex- 
periments on Gd-doped Si in the presence of large mag- 
netic fields show just such a remarkable behavior |Io|].A 
systematic search for these effects in other ferromagnets 
may reveal that they are a universal signature of metallic 
ferromagnetism, which would lend support to the mod- 
els discussed here. The magnitude of the effects how- 
ever, just as the magnitude of the magnetoresistance, is 
likely to vary widely from material to material. As dis- 
cussed in Sect. VII, there should be a definite relation 
between the magnitude of the magnetoresistance and the 
change in optical absorption with magnetic field, if in- 
deed the dominant effect is the lowering of effective mass 
with spin polarization, as suggested by our model, rather 
than a change in the relaxation time. Further discussion 
of these issues and the relation of our model with the 
'double exchange model' [[| is given elsewhere | Q . 

Unfortunately, the model discussed here does not in 
itself contain the physics of the high energy degrees of 
freedom from where the extra optical spectral weight that 
appears at low frequencies gets transfered from. It would 
be of great interest to have a model that would describe 
this high energy physics and give rise to the Hamiltonian 
considered here as an effective Hamiltonian for its low 
energy degrees of freedom. 

Finally, photoemission experiments, especially if angle- 
and spin-resolved, should be able to provide essential 
clues on the validity of the model discussed here. An 
increase in the Fermi velocity of electrons of at least one 
spin orientation as function of increasing magnetic field 
or decreasing temperature (below T c ) would be expected 
to also be a universal feature of metallic ferromagnets if 
the models discussed here are applicable. 

There has recently been substantial interest in re- 
analyzing the problem of metallic ferromagnetism p4] , p5| , 
and there seems to be a consensus that this old problem 
is still not well understood despite the practical successes 
of spin-density- functional theory ||. In particular, recent 
work has suggested that peaks in the density of states are 
the dominant mechanism giving rise to ferromagnetism 
^4j,^5|. It has also been suggested |25| that because fer- 
romagnetism is a strong coupling problem it is not really 
possible to pinpoint the single ultimate cause of metallic 
ferromagnetism, whether peaks in the density of states, 
band degeneracy, or particular electron-electron interac- 
tions. However, the work discussed here and previously 
[plpp]] suggests otherwise. We expect metallic ferromag- 
netism to be always accompanied by the transfer of opti- 
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cal spectral weight discussed above, as well as by a change 
in the quasiparticle dispersion. Establishing experimen- 
tally that these phenomena do not occur would prove 
the invalidity of the model discussed here. Concerning 
band degeneracy, if ferromagnetism is found in a system 
where conduction clearly occurs through non-degenerate 
bands, such as metallic hydrogen [^6), it would establish 
that band degeneracy is not essential and lend further 
support to the model discussed here. 
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FIG. 2. Ground state phase diagrams, a = ji/ja- j = ji- 
P, PF and F denote paramagnetic, partially ferromagnetic 
and fully polarized ferromagnetic regions respectively. In 
the absence of exchange (k — 0), full polarization cannot be 
achieved for band fillings between 1/4 and 3/4; that portion 
of the phase boundary is indicated by a dashed line. As k 
increases the region where full polarization can be achieved 
increases and for k > 0.5 it covers the entire phase diagram 
except for the point n = 1. 

FIG. 3. Magnetization versus temperature for band filling 
n = 0.5, k = 0, jl = 0. The numbers next to the curves give 
the values of ji = j. For j < 1.6 re-entrant ferromagnetism 
occurs in this case. 

FIG. 4. Effective mass ratio, or inverse of bandwidth, ver- 
sus temperature for the parameters of Fig. 3. (a) Majority 
spins; (b) minority spins. The numbers next to the curves 
give the values of ji — j. 



FIG. 5. Effective moment (Eq. (40)) versus temperature 
for the parameters of Fig. 3. 

FIG. 6. (a) Magnetization versus temperature and (b) ef- 
fective moment versus temperature for n — 0.5, k — 0, 
a = Ji/j2 = 0.5. The values of j'2 = j are given next to 
the curves. 

FIG. 7. Effect of exchange splitting on the behavior of (a) 
magnetization and (b) effective moment versus temperature. 
For (a), the curves of increasing steepness correspond to the 
cases (i) - k) = 2, k = 0; (ii) j/(l -k)=2,k = 0.25; (iii) 
- k) = 2, k = 0.5; (iv) j/(l -k) = 4,k = 0.5. 
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FIG. 8. Optical conductivity (arbitrary units) versus fre- 
quency for the parameters of Fig. 3. The full lines give the 
optical conductivities at temperatures T/T c — 0.2, 1 and 3, 
and the dashed lines the corresponding ones in the presence of 
a magnetic field h — H/D — 0.02. The effect of the magnetic 
field is largest at T c , and becomes very small for T/T c = 3. 

FIG. 9. Optical weight n/m* given by the low-frequency 
integral of the conductivity versus temperature in the pres- 
ence of a magnetic field h = H/D (numbers next to the lines), 
(a) Total, (b) for spin up (full lines ) and spin down (dashed 
lines) respectively. The parameters are those of Fig. 3. 

FIG. 10. Magnetoresistance versus temperature for the pa- 
rameters of Fig. 3. 

FIG. 11. Angle-resolved photoemission spectra for the pa- 
rameters of Fig. 3, giving rise to T c / D — 0.081. We assume a 
bandwidth D = 0.3eV, so T c = 280K. Full lines correspond 
to temperature T/T c = 1.2, dashed lines to T/T c = 0.8. The 
numbers next to the full lines give the values of — H (in 
meV). (a) Majority spins; below T c , the Fermi velocity in- 
creases and the peaks disperse faster, (b) Minority spins: the 
Fermi velocity decreses upon ordering and the peaks disperse 
slower, (c) Non-spin-resolved: two peaks appear below T c , 
one disperses faster and one slower than the peaks above T c . 
Similar results would be found at fixed temperature under 
application of a magnetic field. In particular, application of 
a magnetic field h = H/D = 0.03 at T/T c = 1.2 gives results 
almost identical to those shown in the figure for T/T c — 0.8. 
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